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Abstract. - We analyze the effect of a quantum error correcting code on the entanglement 
of encoded logical qubits in the presence of a dephasing interaction with a correlated environ- 
ment. Such correlated reservoir introduces entanglement between physical qubits. We show 
that for short times the quantum error correction interprets such entanglement as errors and 
suppresses it. However for longer time, although quantum error correction is no longer able to 
correct errors, it enhances the rate of entanglement production due to the interaction with the 
environment. 



Quantum error correction (QEC) [1,2] has been introduced to perform quantum infor- 
mation processing [3] in the presence of noise due to the interaction between quantum bits 
and environment. The basic idea is to encode each logical qubit in a redundant way on a 
set of physical qubits and to periodically acquire information on the errors that affected the 
system but not on the quantum state of the system itself. Such techniques have been devel- 
oped to deal with independent errors on individual physical qubits due to the interaction of 
each physical qubit with its own reservoir. In several physical situations however the presence 
of correlated reservoirs [4] can result in non-trivial effects. For example it has been shown 
that the interaction of two subsystems with a finite temperature common bath of harmonic 
oscillators can, for short times, induce entanglement between the two subsystems initially in 
a product state. This is possible when the environment has some spatial correlations [5], as 
often occurs in solid state physics, leading to an effective interaction between the two subsys- 
tems. This effect is also present for noisy baths in the Markovian regime [6]. The dynamics 
of the entanglement rate in the presence of decoherence was also studied [7] . 

Quantum error correction has been analyzed in the presence of correlated environments in 
Ref. [8] . In the present work we address the issue of the effects of QEC on the entanglement 
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between logical qubits. It is reasonable to expect that the entanglement induced by the corre- 
lated bath between physical qubits will modify the encoded state in a way that is interpreted 
by the QEC procedure as error and therefore corrected. However, when such entanglement 
becomes sufficiently large the protocol may be not able to correct it. It is therefore interest- 
ing to study how entanglement is modified by the application of QEC. In the following we 
will show that, although QEC is unable to correct such errors, it can enhance the generation 
of entanglement in a pair of logical qubits with respect to the entanglement induced by the 
environment on a pair of physical qubits. 

The model we consider in this work, the same as in [4], consists of a register of quantum 
bits interacting with a common environment, modelled as a bath of harmonic oscillators. The 
bath - qubit interaction is described by the following Hamiltonian 



(1) 



where £j(t) — w [Aj lTO (w)a TOia ,+/ic]. In the previous expression A Jim (w) denote the coupling 
constants between the jth qubit and the oscillator at frequency w in the mth bath with 
corresponding annihilation (creation) operator a m ^(al n u ). 

In the following we will concentrate our interest on the register dynamics i.e. on the reduced 
density operator p(t) = Tr E[U(t)p(0)U^(t)] where Tr^ denotes the partial trace performed on 
the environment degrees of freedom. The resulting density matrix can be written in a compact 
form as [8] 



p(t) = $[U r p(O)0] 
where $ is a map (a super-operator) defined as 



(2) 



= exp 



3,3' 



(3) 



In the above expression we used the convention that a bar over an operator means that it 
acts on the density matrix from the left and r,y(i) = {(j}j{t)4>ji (t)) , with <pj{t) — J £j(t')dt', 
involves correlations of the bath at different times. Finally the unitary evolution U r is given 

by 



U r = exp 



(4) 



where the quantity Vjji = 2Rc ~^2 m u X m ,j ( w )^* n y ( w ) /w is non zero only if the same reservoir 
is coupled to different qubits. As we will see in the following, the unitary operator U r is re- 
sponsible for the creation of entanglement, while the super-operator $ describes the dephasing 
of the off diagonal elements and is present also without spatial correlations of the environment. 

As a measure of entanglement between two qubits we will use the concurrence defined by 
Wootters [9] . If p is the density matrix of the global system of two qubits, let us define p = 
g v ®a v p* u v ®a y and R = pp. The concurrence is then defined as C — max{0, Ai — A2 — A3 — A4} 
where Ai are the square roots of the eigenvalues of R labelled in decreasing order. Using this 
definition we obtain for Bell states 1^) = (|00) ± |H))/\/2 the result C(t) = e - 4 ( r n+ 2r i2+r22) 
while for |V' ± ) = (|01) ± |10>)/V2 the concurrence decays as C(t) = e - 4(ril - 2Tl2+r22 \ It is 
easy to see that if the environment has special symmetries, for instance if is equal for 
different i and j then the subspace spanned by |i/' ± ) is decoherence free [4,10]. Note that the 
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typical decoherence times are of the order of dt/dTi. If we choose \ij) x , eigenstates of a x as 
initial states, it is possible to see that the qubits become entangled if V12 > dTi/dt. In fact 
the entanglement oscillates with a frequency AV12 and is damped with a rate proportional to 
dTi I dt (see figure . 




Fig. 1 - Time evolution of fidelity and concurrence between physical qubits, averaged over initial 
states, in the absence of QEC. The time scale is chosen so that V = 1 and F(t) ~ O.lt 

Let us now introduce a QEC protocol. In the following we will encode each logical qubit 
into three physical qubits as follows |5) — * (000)^ and l) — > 1111)^. where 10}^ and ll)^ are 
eigenstates of a x . This code will protect the logical qubits against one single phase error on 
individual qubits [1]. When VT, T(T) <§C 1, where T is the time interval between two quantum 
error corrections, following the same approach of [8] we can write the following master equation 
for the continuous evolution of the density matrix of the two logical qubits 

-£ = -(7i + li)p + Tioipoi + 72^P^ ( 5 ) 

where 

7i = I £ (T 2 V?f + 2T l' ( T ) + T n (T)T rr (T)) . (6) 

]>j' 

as in Ref. [8]. In the above equation we have to take j = 1, 2, 3 for i = 1 and j — 4, 5, 6 for 
i = 2. The rates % <§C dTij/dt are the same found in [8] for the single qubit case. 

Let us now turn our attention to the time evolution of the entanglement between logical 
qubits in the presence of QEC in the regime in which the above master equation holds. As a 
first example we consider the initial state |00). In this case the only non- vanishing elements 
at time t are those in the main diagonal, and they decay exponentially with decay constants 
equal to 71,72 and 71 +72. If instead we start with a Bell state \<j> + ) of logical qubits the 
non-zero elements at time t are those on the two diagonals. Such state evolves towards an 
incoherent superposition of two Bell states whose concurrence decays as C(t) = e -2 ( 7l+72 )* 
and the fidelity with respect to the initial state is F(t) = h(l + e -2 ^^ 72 '*), where, for an 
initially pure state \ip(t = 0)) the fidelity is defined as F(t) = (^f(0)|p(i)|^>(0)). From these 
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results it is evident that the fidelity and the concurrence, in the presence of QEC, decay with 
the same rate. It is interesting also to see what happens to the state . In this case the 
density operator does not evolve: pit) = p(Q). This does not follow from the decoherence 
free subspace phenomenon but rather from the fact that QEC freezes the unitary evolution by 
means of repeated measurement, as in the quantum Zeno effect [11], and conditional dynamics 
depending on the measurement outcome. 

In order to obtain information about the global properties of the map we will consider 
states with fixed initial concurrence, which are of the form: 

|*) =costf |0„i0„ 2 ) +sintf|l n ii n2 ) (7) 

with 

|6 m ) = cos^|0)+sin^e^ |l) (8) 

|l ra ) = cos||l>-sm^e-^|6) (9) 

The concurrence of |^) is sin 2$. In the following analysis we will average over -d 2 , ¥>i> ¥?2- 
We can see that averaging over all initial product states $ = 0, ir/2 one finds that the con- 
currence is always zero and the fidelity decays with a rate proportional to 7^. For partially 
or maximally entangled states the entanglement decays monotonically again with decay con- 
stants 71+72. This means that QEC suppresses the effective interaction between logical qubits 
due to the presence of a correlated environment. In other words the bath does not induce en- 
tanglement between logical qubits and the entanglement initially present decays. This implies 
that for small T there is not creation of entanglement since the QEC protocol destroys all the 
correlations between physical qubits of different logical qubits. The reason for this can be seen 
in a qualitative way: for a state |00) = 1000)^ 1000}^ the operator U r = 1 +iJ2ij Vijt cr l a i U P 
to first order creates superpositions like 1000)^ |000) x + IllO)^ 1000)^ + • • • which is entangled. 
After the corrections one gets a mixture of |000) x |000) x and llll)^ |000) x which is no longer 
entangled. Not surprisingly then QEC inhibits the production of entanglement as this is seen 
by the protocol as an error. 

This may no longer be true if the time T is comparable to the period of oscillation of 
entanglement without QEC. In this case the environment may have time to create enough 
entanglement to be interpreted as a property of the initial state and amplified by QEC. If 
the time T between two corrections is not short the approximations which lead to the master 
equation JJjJ are no longer valid and we must use the exact map which links p(t) and pit + T) 
after the free time evolution and the QEC. In terms of Kraus operators such map can be 
written as 

p c (t + T) = J2M k p(t)Ml (10) 

k 

where £ fc M fe f M fe = 1 . In our case we found the following seven Kraus operators: 

Mo — mol 

Ml = TOlCTzl 

M 2 = m,2cr X 2 

M 3 = m 3 a x ia x2 

M4, = mi(a x i + <r x2 ) 

M 5 = to 5 (]1 +<7 x ia x2 ) 

M 6 = \m 6 \(l ±i<j xl a x2 ) 



(11) 
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Such Kraus operators are linked to the (T) rates, however the expressions which explicit 
such dependence are lengthy and not of immediate reading. A more straightforward physical 
picture of the action of the map is instead gained by considering the plot of the coefficients 
mi as functions of T (see FigEJ. We supposed that Tij(T) = T(T) ~ 0.1T from which it 
follows that mi = m-z. Notice that in general 7714 = TO5 and that the sign in Mq is plus 
or minus depending on the sign of Im(me). Figure shows clearly various regimes for the 




Fig. 2 - Relative weights of the Kraus operators, i.e. the quantities mj, as functions of T. For 
simplicity we assume Tij(T) — T(T) ~ 0.1T. From this it follows that mi = 7712. 

evolution of the logical qubits. Indeed for VT <C 1 only too and mi are non zero, leading to 
the dissipative master equation (J5J. As we have seen in this regime there is no creation of 
entanglement because Mi^ are single qubit operators. When T is close to the value ir/8V 
there are also contributions from the two-qubit operators M4 5,6 but the dissipative effect of 
-^1,2,3 is predominant and again we verified that the entanglement production is zero. The 
interesting regime is for T ~ when the main contribution comes from M 6 . As expected 
in this case there is creation of entanglement that oscillates between zero and a maximum 
value (which will gradually decay). Indeed for the particular choice T = (see figure 01 the 
created entanglement is maximum and the QEC amplifies it. The other non zero dissipative 
contribution is Mi and it is present only if Ti is non zero. On the other hand, if Ti is negligible, 
Mq is the only Kraus operator, and therefore the evolution is unitary even in the presence 
of QEC. Under this condition this evolution induces the transformation |6) — * |l) and vice 
versa. For larger T one again finds regimes with no entanglement production. 

When the map I|1U[) is applied (for I\ = and T = yr?) to the state 00) then one sees 
that it makes a sort of periodic oscillation: 

1 06) -► 1 66) + i 1 11) -> |ii) -► 1 66) — i |ii> -► 1 66) (12) 

Of course if T is not zero then the other Kraus operators also contribute and so this oscillation is 
damped. This type of oscillation is not restricted to initial product states but also to partially 
entangled states. As a consequence it is possible to create entanglement that can be distilled 
to obtain maximally entangled states or even extracted from the system. 
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Fig. 3 - Mean fidelity and concurrence between logical qubits in the presence of QEC as functions of 
the number of applications of the map, for T = 7r/4V, V = 1 and T(t) ~ O.li, the same parameters 
used in figure Q 



Figure |21 shows that the maximum average concurrence achievable is around 0.4. It is 
important to underline that this is an average value: there are states, like |00) , that do not 
evolve and so there is no production of entanglement while there are states for which the cre- 
ated entanglement is more than the average value. This is the case of initial product states of 
the computational basis for which the production of entanglement is maximum. For example, 
let us compare the initial states that lead to the maximum production of entanglement with 
and without QEC. In the presence of QEC the best case is |00): the entanglement is 0.92 
after just one application of QEC. On the other hand the maximum entanglement reached 
in the case without QEC is only 0.74 for the state |00) . These results are insensitive to the 
change — > 1. 

In summary, we studied the effect of QEC on the entanglement between logical qubits in 
the presence of correlated noise. We found that when the time interval T between consecutive 
applications of QEC is small QEC is helpful in preserving the state of the qubits and indeed 
entanglement decays with a smaller rate. As expected, in the absence of QEC there is pro- 
duction of entanglement between physical qubits, due to the correlations introduced by the 
environment, while in the presence of QEC this phenomenon is inhibited. The interesting new 
effect we have discovered is that when T is not small then entanglement may be generated 
with a bigger rate than in absence of QEC. In this scenario we are in a situation in which the 
entanglement induced by the unitary dynamics generated by the correlated bath is enhanced 
by the QEC protocol which prevents the dephasing effects of the coupling with the bath. Note 
that such enhancement in the production rate of entanglement is achieved by means of local 
measurements and conditional local unitary operations on the logical qubits, in other words 
the entanglement is not induced by joint measurements on the pair of logical qubits. 

Our results suggest a new strategy to enhance the rate of entanglement production for 
interacting qubits in the presence of decoherence in a more general scenario. Consider a system 
of interacting subsystems in the presence of a decohering dynamics. Such situation might be 
due either to a coupling of the subsystem to a common reservoir, like in our case considered 
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above, or to a direct mutual coupling between subsystems in the presence of independent 
environments. The entanglement induced between subsystems by the unitary dynamics is 
destroyed by the dissipative one. By enlarging the Hilbert space, in our case by increasing the 
number of qubits, it is possible to identify separate larger subsystems such that by means of 
local measurements and local conditional unitary operations on such subsystems the amount 
of entanglement due to the mutual interaction is increased with respect to the entanglement 
it would have been generated between the original subsystems. In other words we conjecture 
that the strategy outlined above is not restricted to QEC but it is a more general technique. 
The issue is, given an effective interaction in the presence of decoherence, to find optimal 
subsystems and projection protocols in order to maximize the production of entanglement. 

In some sense our protocol can be described as a generalization of the one proposed in Ref. 
[12]. In [12] the entanglement production is optimized, in the presence of a direct interaction 
and in the absence of decoherence, by means of local operations and ancillas. In our QEC 
protocol we use some sort of ancillary system to enlarge the Hilbert space, although in this 
case there is not a sharp distinction between qubits and ancillas. Furthermore, we make use 
also of local projections on the enlarged subsystems and conditional dynamics. 

We gratefully acknowledge many helpful discussions with D. Averin. This work was sup- 
ported by the EU (IST-SQUBIT, HPRN-CT-2002-00144; QUPRODIS: IST-2002-38877). 
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